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We consider a simple version of a cyclic adiabatic inver- 
sion (CAI) technique in magnetic resonance force microscopy 
(MRFM). We study the problem: What component of the 
spin is measured in the CAI MRFM? We show that the non- 
destructive detection of the cantilever vibrations provides a 
measurement of the spin component along the effective mag- 
netic field. This result is based on numerical simulations of 
the Hamiltonian dynamics (the Schrodinger equation) and the 
numerical solution of the master equation. 
PACS: 03.65.Ta, 03.67.Lx, 76.60.-k 



I. INTRODUCTION 

Magnetic resonance force microscopy (MRFM) is ap- 
proaching its ultimate goal: a single spin detection [1-3]. 
The most promising approach to single spin detection is, 
probably, cyclic adiabatic inversion (CAI) [1] . In this ap- 
proach, the magnetic moment of the sample changes its 
direction adiabatically following the effective magnetic 
field. The CAI of the spin may act as an "external force" 
driving the resonant vibrations of the cantilever or it may 
affect the frequency of the cantilever vibrations driven 
by another source (e.g. the modern "OSCAR" technique 

The fundamental question which arises in MRFM 
single-spin measurement is the following: What com- 
ponent of the spin is measured by this technique? In- 
deed, in a simple geometry the cantilever tip oscillating 
along the z-axis interacts with the ^-component of the 
spin and, consequently, is expected to measure the spin 
^-component. From the other side, adiabatic inversion 
assumes that the approximate integral of motion is the 
spin component along the effective magnetic field which 
rotates in the x — z plane. Thus, one might expect that 
the cantilever measures the spin component along the 
effective magnetic field in the rotating reference frame. 

In this work we consider the macroscopic cantilever 
as itself the measuring device interacting with an envi- 
ronment. We assume that the influence of an additional 
(e.g. optical) device which detects the cantilever vibra- 
tions is small. This corresponds to the current MRFM 
technique. In Section II we discuss the quantum dynam- 
ics of the quasi-classical cantilever which describes a gen- 
eration of the Schrodinger cat state associated with two 



possible projections of the spin. In Section III we include 
an interaction of the cantilever with an environment in- 
herent to any measurement processes. The latter leads 
to the decoherence of the two possible trajectories of the 
cantilever due to interaction with the environment. 



II. HAMILTONIAN DYNAMICS 

We considering the simple setup shown in Fig. 1. 

The ferromagnetic particle with a magnetic moment m 
is mounted on the cantilever tip. The permanent mag- 
netic field, B , points in the positive z-direction. A ro- 
tating rf field in the x — y plane, B\ ~ exp[i(ujt — ip(t))], 
is resonant with the spin precession around the z-axis. 
The frequency modulation of B\ causes the CAI of the 
spin. Under resonant conditions, when the period of the 
cantilever vibrations matches the period of the CAI, the 
amplitude of the cantilever vibrations is expected to in- 
crease providing the detection of the spin. 

The quantum Hamiltonian of the system in the rotat- 
ing frame ( in terms of dimensionless parameters) can be 
written as 



H={ P 2 Z + z 2 )/2 + <p(r)S z - eS x - 2 V zS z 



Here 



Pz 



P z /P q , z = Z/Z q , 



(1) 



(2) 



S is the electron spin operator, e = ^Bi/lJc, ip — dip/ 1 dr, 
7] = gF/2F q , P z and Z are the operators of the ef- 
fective momentum and coordinate of the cantilever tip; 
7 = gfiB/h is the spin gyromagnetic ratio (absolute 
value); lo c is the cantilever frequency; F is the mag- 
netic force between the ferromagnetic particle and the 
spin when the cantilever tip is at the origin (z = 0). The 
origin is chosen at the equilibrium position of the can- 
tilever with no spin; r = uj c t is a dimensionless time. 
The units of the coordinate, momentum, and force are 
given by 



Z q = (hiv c /k c ) 



1/2 



P q = h/Z q , F q = k c Z q 



(3) 



where k c is the cantilever spring constant. Note, that we 
treat an electron spin of a paramagnetic atom whose di- 
rection is opposite to the direction of the atomic magnetic 
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moment. We assume in (1) that the transverse magnetic 
field points in the negative x-direction of the rotating 
frame. 

With respect to actual "reading" devices we consider 
a realistic scenario for the MRFM technique which in- 
volves non-destructive measurements of the amplitude, 
frequency, and phase of the cantilever vibrations, for ex- 
ample by using a fiber-optic interferometer operating in 
the infrared region. We assume that the optical detection 
of cantilever vibrations does not influence significantly 
the cantilever-spin dynamics. (In practice, this means 
that the disturbance caused by the optical radiation is 
smaller than the thermal noise of the cantilever.) 

In this section we do not consider the interaction with 
the environment which provides the measurement itself. 
(See also [4].) Thus, we use the Schrodinger equation, 



(4) 



for computer simulations of the cantilever-spin dynam- 
ics. In the z — S' s -rcpresentation, the wave function, \&, 
is a spinor. It contains two components, &(z, 1/2, r) and 
4»(z, — 1/2, r), which correspond to the two possible val- 
ues of S z . Using the expansion over the eigenfunctions, 
u„, of the oscillator Hamiltonian, [p z + z 2 )/2, we write 
these two components of the cantilever-spin wave func- 
tion in the form, 



*(*,l/2,r) = VA,(rK, *(z,-1/2,t) = VB n (rK, 



n=0 



n=0 



(5) 

and derive equations for the amplitudes, A n and B ni 
iA n = (n+l/2+ l f/2)A n -(r ] /V2)(V^A n - 1 +v / n + lA n+1 ) 



-(e/2)Bn, 



(6) 



iB n = {n-l/2+£/2)B n +(r ] /V2)(V^B n - 1 +V^TTB n+1 ) 
-(e/2)A n . 

The initial conditions describe the quasi-classical state of 
the cantilever tip and a spin which points in the positive 
^-direction, 



The value 77 = 0.3 can be achieved in current MRFM 
experiments [1-3]. The parameters for the transverse 
magnetic field have been chosen to satisfy two conditions: 
i) the condition of the CAI, \(p\ -C e 2 , and ii) the effec- 
tive magnetic field produced by the cantilever vibrations 
on the spin is small in comparison with the amplitude of 
the rf field: 27/1 <C e. We consider the results of the 
computer simulations reliable if they do not change with 
an increase in the number of basic functions, u n . 

To describe the quasi-classical cantilever, we took the 
initial average energy (E(0)) = \a\ 2 3> 1. The number 
of basic functions, u n , needed to provide reliable results, 
increases with the average energy. So, we can not take \a\ 
too big. As we study the driven oscillations of the can- 
tilever, our results did not show a significant dependence 
on the initial conditions. 

The main results of our simulations are the following. 
The wave function of the cantilever-spin system, which 
is initially a product of the cantilever and spin parts, 
quickly becomes entangled. The probability distribution 
to find the cantilever at the point z at time t, 



P(z,t) = |*(z,l/2,r)| 2 + |*(z,-l/2,r)| 2 



(9) 



splits into two peaks, "big" and "small" peaks. (Sec Fig. 
2.) When the peaks are separated, the wave function of 
the cantilever-spin system can be represented as a sum 
of two spinors, 

*(z,s,t) = * (1) (z,s,t) + * (2) (z,s,t), (10) 

where the upper indices "1" and "2" refer to the big 
and the small peaks, correspondingly. It was found with 
the accuracy to 1% that both spinor wave functions, 
$W(z,s,t) (k = 1,2), can be represented as a product 
of the cantilever and spin functions, 



HfW{z,8,t) = rW{z,t)x W {8,t), 



(11) 



where ( s j T ) describes the spin which points in the 
direction of the external effective field, (e, 0, — <p(r)), and 
X^( s j r ) describes the spin which points in the opposite 
direction. The ratio of the probabilities for the big and 
the small peak is determined by the initial angle between 
the external effective magnetic field and the spin, 



A„(0) = K/V^!)cxp(-| a | 2 /2), B„(0) = 0, (7) J ' \R {2 \z,r)\ 2 dz J \R^(z, r)\ 2 dz = tan 2 (6/2), 



a=[(z(0))+i( Pz (0))}/V2. 

In our computer simulations we used the following pa- 
rameter values, 

77 = 0.3, e = 400, 93 = -6000 + 300r if t < 20, (8) 



and ip = 1000sin(r - 20), if t > 20. 



(12) 

where 8 is the initial direction of the external effective 
magnetic field (tan© = —e/ip(0) = 1/15). If the ini- 
tial conditions describe a spin which points, for example, 
in the positive x-direction (A n (0) = B n (0j), our simula- 
tions reveal two peaks with approximately equal ampli- 
tudes. Thus, the Hamiltonian dynamics clearly indicates 
that the quasi-classical cantilever will measure the spin 
component along the effective magnetic field. Certainly, 
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in the frames of the Hamiltonian approach we cannot 
describe the measurement itself: the coherence between 
the two cantilever peaks does not disappear. In other 
words, the Schrodinger equation describes the macro- 
scopic Schrodinger cat state of the cantilever without ef- 
fects of dccohcrcncc. 



Next, we solve numerically the system of equations for 
the amplitudes, A^ s m {T), 



(15) 



[itp(r) («' - s) +(3/2 - (n+m+ 1)D0 - i(n - ro)] A^ n (r)- 



III. MASTER EQUATION 

In the previous section, we have presented indications 
that the cantilever "measures" the spin component along 
the direction of the effective magnetic field. In this sec- 
tion we describe the measurement process. During the 
measurement process the coherence between two can- 
tilever trajectories disappears. It means that the reduced 
density-matrix of the cantilever-spin system becomes a 
statistical mixture representing two possible trajectories 
of the system. The main question we are going to answer 
is the following: Does the cantilever, which interacts with 
the environment, measure the spin component along the 
effective magnetic field? 

To answer this question, we studied the dynamics of 
the cantilever-spin system using the master equation. 
Our purpose is not just to simulate the expected exper- 
iment but rather to present a qualitative verification of 
the conclusion obtained in the previous section. That is 
why we consider the simplest "ohmic" model of the envi- 
ronment in the high-temperature approximation [5]. In 
this approximation the environment is described as an 
ensemble of harmonic oscillators. The number of oscilla- 
tors per unit frequency is proportional to the frequency 
in the region below the chosen "cutoff" frequency, f2, and 
ksT 3> Kl. The master equation for the density matrix, 
/?, in the high-temperature approximation, is 



dr 




dz' 



~- 2 ^ 2 ~n as) 



iijs'V^A^^t) - tr ]S W2m + 2A^ m+1 (T) + 
i V sV2HA s / hm (T) + ir]sy/2n + 2A s 4 hm (T)+ 



D/3y/m(n + ^A^^+D^^m + l)<l 1>m+1 (r)+ 



(D + l/2)/V(" + l)(m+lK'; M „ +1 (T) + 



(D-l/2)M^li, m -iM- 



(D-l/2)^^(^T)A^ m (r)- 



(D + l/2)^(n + l){n + 2)A s n ^ m {r) 



(D l/2)^m( m -l)A^ m _ 2 (r)- 



(D + l/2)|V(m + 2)(m+l)<f ro+2 (T)- 



i^ m s '(r)-A-^'(r)}. 



~2 {Z - 



-2irj(z's' — zs) + i(p(s' — s) 



Pss'(z, 



Z ,T 



e 

- l 2 



Pss'(z,z',t) - p Ss '(z,z',T) 



Here s, s' = ±1/2, s = — s, s 1 = —s', D = kBT/hu> c , 
[3 = l/Q, where Q is the quality factor of the cantilever. 
Again, we use the expansion over the eigenfunctions, u ni 

PS , S ,(Z,Z',T) = Y J A ninij)u n {z)u* m {z'). (14) 



Below we describe the results of our computer sim- 
ulations for the values of parameters in (8). First, 
setting (3 = D = we obtain the density matrix, 
p s ^s' (z, z', t), which exactly corresponds to the wave func- 
tion, ^(z, s,t), derived from the Schrodinger equation. 

The initial density matrix is represented as a product 
of the cantilever and spin parts, 

Ps , s/ (z,z',0) = *(z,l/2,0)**(z', 1/2,0) (J (16) 

The wave function, ^(z, 1/2,0), describes the quasi- 
classical state of the cantilever, 



f(z,0) = ^A„(0K(z). 



(17) 



n=0 
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The values A n (0) are given in (7). The initial values, 
An m(0)> m (15) can be easily found from (16). 

For t > 0, the density matrix describes the entangled 
state which cannot be represented as a product of the 
cantilever and spin parts. The initial peak of p SiS > (z,z', r) 
splits into two peaks which are centered along the di- 
agonal z = z' , and two peaks centered at z ^ z' , off 
the diagonal. The density matrix can be represented ap- 
proximately as a sum of four terms corresponding to four 
peaks, 

p s ,Az,z',r) = P^+P%+P%+P%, (18) 

where we omit variables, z,z',t. The matrices, and 
p( 2 ' , describe the "big" and "small" diagonal peaks; p^ 
and p^\ describe the peaks centered at z ^ z' . 
As an illustration, we show in Fig. 3 the quantity, 

IPl/2,1/2^, z', t) + p_ 1/2 _ 1/2 (z, z', t)|. 

We have found that with accuracy to 1% the density 
matrix, p^\,{z, z', r), can be represented as a product of 
the coordinate and spin parts, 

P^,(z,z',t) = R^(z,z',t)x^,(t), (19) 

where X^ 1 ]/ (t) describes the spin which points in 
the direction of the external effective magnetic field 
(e, 0, — <p(t)). A similar expression is valid for 
p^\, (z,z',t). But in this case, X^l>( T ) describes a spin 
which points in the opposite direction. 

First we note that in order to describe the measure- 
ment process (the decoherence) , we have to consider an 
ensemble of quasi-classical cantilevers with the same ini- 
tial conditions. At the same time, we are considering 
driven oscillations of the cantilever. So, the result of our 
simulations qualitatively does not depend on the initial 
conditions of the cantilever. Second, as we already men- 
tioned, we are going to verify qualitatively the conclusion 
derived in the previous section rather than simulate the 
expected experiment. Thus, we choose the values of pa- 
rameters which help us to save a computational time. 
Namely, we choose a relatively small (but still quasi- 
classical) values for the initial energy of the cantilever, 
and a relatively small value for the thermal parameter 
D (without violating the high-temperature approxima- 
tion which requires D 3> 1). The small initial energy 
of the cantilever allows us to reduce a number of basis 
functions, u n (z). A relatively small value of D allows us 
to observe four well separated peaks at relatively small 
values of time, r. 

The initial uncertainty of the cantilever position is, 
5z = l/\/2. Due to thermal diffusion, the uncertainty 
of the cantilever position increases with time. Thus, we 
have two effects: i) the increase of the amplitude of the 
driven cantilever vibrations (similar to the Hamiltonian 



dynamics) and ii) the increase of the uncertainty of the 
cantilever position due to thermal diffusion. If the sec- 
ond effect dominates, the two positions of the diagonal 
peaks (i.e. peaks centered on the line z — z') become 
indistinguishable. In this case, one cannot provide a spin 
measurement with two possible outcomes. 

We have found that peaks centered on the diagonal 
retain the main properties described by the Hamiltonian 
dynamics. The density matrix, p^\, (z,z',t), for k = 1,2 
can be approximately represented as a product of the 
cantilever and spin parts. The spin part of the matrix 
describes the spin which points in the direction of the 
external effective magnetic field (k = 1) or in the opposite 
direction (k = 2). 

Next, we discuss the two peaks centered at z ^ z' . As 
an illustration, Figs. 4 and 5 show the contours of the 
quantities, 

\Pl/2.1/2(z, z', t) + p_i/2,-l/2(z, z', T)|, 

and 

|Pi/2 -1/2(2, z', t) + ,0-1/2,1/2(2, z, r)|, 

given in logarithmic scale. One can see the peaks cen- 
tered at z 7^ z' as well as at z = z'. The peaks centered 
at z 7^ z' describe the coherence between the two can- 
tilever positions. The amplitude of these peaks quickly 
decreases due to the decoherence. Thus, the master equa- 
tion explicitly describes the process of measurement. The 
coherence between two cantilever trajectories (the macro- 
scopic Schrodinger cat state) quickly disappears. As a 
result, the cantilever will "choose" one of two possible 
trajectories. Correspondingly, (depending on the can- 
tilever trajectory) the spin will point in the direction of 
the effective magnetic field or in the opposite direction. 

IV. CONCLUSION 

We have studied the quantum dynamics of the 
cantilever-spin system in a simple version of the cyclic 
adiabatic inversion (CAI) magnetic resonance force mi- 
croscopy (MRFM) . In this version, the spin experiences a 
CAI under the action of the external phase modulated rf 
magnetic field. If the frequency of CAI matches the can- 
tilever frequency, the amplitude of the cantilever vibra- 
tions increases allowing single-spin detection. We have 
studied the problem: Which component of the spin is 
measured by the cantilever? We argue that one will mea- 
sure the component of the spin along the direction of the 
effective magnetic field providing non-destructive detec- 
tion of the cantilever vibrations. This result was first 
derived using computer simulations of the Hamiltonian 
dynamics (the Schrodinger equation). Then, it was con- 
firmed by the numerical solution of the master equation. 
We have considered the case when the amplitude of the 
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driven cantilever vibrations was greater than the thermal 
noise. In this case, the phase of the driven vibrations 
depends on the spin component along the direction of 
the external effective magnetic field. Thus, detecting the 
phase of the cantilever vibrations one can measure the 
spin component along the effective magnetic field. 

We should mention that the direct relation between 
the cantilever trajectory and the direction of the spin has 
been verified for a transient process in the CAI MRFM. 
Our computer capabilities do not allow us to check this 
relation for the stationary cantilever vibrations at r 3> Q- 
Also, we completely ignored the direct interaction be- 
tween the spin and the environment. We are now inves- 
tigating this interaction. 
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FIG. 1. MRFM setup. 



FIG. 2. The probability distribution, P(z,t), for the can- 
tilever position, in the logarithmic scale. The values of pa- 
rameters are: e = 400 and r\ — 0.3. The initial condi- 
tions are: (z(0)) = —20, (p z {0)} = (which corresponds to 
a = -10-^2). 



FIG. 3. Three-dimensional 

plot of log |pi/ 2 ,i/ 2 (z,z',r) + /o_i/ 2 ,_i/2(z,z',t)|, in the log- 
arithmic scale. The values of parameters are: e — 400, 
r\ = 0.3, (3 = D = 0. The initial conditions are: (z(0)) = —4, 
(pz(0)) = 0. The colors in Figs. 3-5 correspond to the fol- 
lowing values for the logarithm of the density matrix: white 
(<-16), red (-16,-12), green (-12,-8), blue (-8,-4), yellow (>-4). 



FIG. 4. The 

contours for log |pi/ 2 , 1/2(2, 2', t) + p_ 1/2 ,-1/2(2, 2', r)|. The 
values of parameters are: e = 400, (3 = 0.001, and D — 10. 
The initial conditions are: (z(0)} = —4, (p z (0)) = 0. 



FIG. 5. The same as in Fig. 4, but for 

log \pi/2, -1/2(2, 2', r) +P-1/2, 1/2(2, z',t)\. 
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